Shot Noise in Ballistic Graphene 
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We have investigated shot noise in graphene field effect devices in the temperature range of 4.2-30 
K at low frequency (/ = 600-850 MHz). We find that for our graphene samples with large width 
over length ratio W/L, the Fano factor 5 reaches a maximum 5 ~ 1/3 at the Dirac point and that it 
decreases strongly with increasing charge density. For smaller W/L, the Fano factor at Dirac point 
is significantly lower. Our results are in good agreement with the theory describing that transport 
at the Dirac point in clean graphene arises from evanescent electronic states. 

PACS numbers: 



03 



c^ 



O 

> 

O 



o 



^ 



Recently, graphene, a two-dimensional carbon crystal 
arranged in honeycomb lattice, has attracted tremendous 
attention in the scientific community [1]. Its fundamen- 
tal physical properties open up possibilities jor realizing 
high-speed, ballistic field effect transistors [2] as well as 
for potential applications like spin control devices [3] or 
gas sensors \^. Owing to its unique band structure, low 
energy conduction in graphene occurs via massless Dirac 
fermion quasiparticles leading to very rich and sometimes 
counterintuitive behavior [1]. For example, despite the 
vanishing density of states at the Dirac point where the 
conduction and the valence band touch, the conductiv- 
ity of graphene remains finite. Indeed, it has been pro- 
posed that transport at the Dirac point can be explained 
as propagation of charge carriers via evanescent waves 
[5|. Universal minimum conductivity of ^ has been 
observed around the Dirac point for samples with large 
width over length ratio (W/L) [2|, yielding support to 
this evanescent- wave theory [5||. In addition to the min- 
imum conductivity, current fiuctuations are expected to 
exhibit non-trivial characteristics. A universal maximum 
value of 1/3 for the Fano factor 5^, decreasing with in- 
creasing charge carrier density, is expected for large W/L 
in graphene. However, no experimental results have been 
reported yet [6]. In this Letter, we present a study of shot 
noise in short graphene strips at frequencies / = 600—850 
MHz, which is well above the 1// corner frequency, even 
at the large currents when working far away from the 
Dirac point M. Using a high resolution noise measure- 
ment set-up [8|, we measure shot noise of two-terminal 
graphene devices both as a function of the charge density 
as well as bias voltage. Our results indicate that ballistic 
transport occurs at the Dirac point via evanescent waves. 

The role of disorder, interactions or carrier statistics in 
quantum transport can be probed by shot noise in meso- 
scopic devices [9]. Arising from the granular nature of 
electron charges, shot noise reveals information on funda- 
mental conduction properties of low-dimensional systems 
in a complementary manner to conventional dc transport 
measurements. Shot noise can be detected only when the 
drift time is shorter than electron-phonon energy relax- 
ation time [10, Uj|. Naively, since electrical transport 
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FIG. 1: (a) Experimental set-up for detecting shot noise at 
T = 4.2-30 K. (b) Schematic of the principle of our measure- 
ments in terms of the noise power reflection |r|^. (c) Illus- 
tration of a typical graphene sample fabricated for our shot 
noise study. 



in graphene at submicron scale is believed to be ballis- 
tic, one would expect shot noise to be completely absent, 
Le. the Fano factor d = En=i ^n(l - T^)/ EnLi ^n = 
0, where T^ are the transmission coefficients of the con- 
ducting channels in the scattering matrix formalism [9| 
(when ^ = 0, all eigen values of the transmission matrix 
are equal to 1). However, Tworzydlo et al. [5] have calcu- 
lated that perfect graphene, contacted by highly doped 
graphene leads, should exhibit a Fano factor of 1/3 at 
the Dirac point for high W/L (note that Katsnelson re- 
ported similar sub-poissonian shot noise in [12|). Similar 
calculations using metallic square lattice contact lead to 
the same conclusions [l3|. In Ref. |5|, the authors show 
that the transmission coefficients T^ = cosh~^[7rn-^], i.e. 
they display the same distribution as in coherent disor- 
dered conductors. No clear-cut explanation for this in- 
triguing coincidence between transmission eigenvalue dis- 
tribution for the propagating modes in diffusive systems 
and evanescent states in graphene has been given [5]. 

Our experimental setup is illustrated in Fig. 1(a). In 
our work, we use the sensitive lock-in detection technique 



described in Refs. |8|, [iJ. To enhance the sensitivity of 
noise detection, we use a low-frequency sine- wave modu- 
lation of the current, / = /c^c + ^^sin(c(;t) where Idc ^ ^^• 
In order to avoid external spurious signals, the set-up is 
placed in a Faraday cage and the signal is band limited 
to the range of 600-850 MHz, thus avoiding EMI from 
mobile phone frequencies. We measure the shot noise of 
a tunnel junction (g^ = 1) for calibration. A microwave 
switch is used to alternatively measure the noise from 
the graphene sample and the tunnel junction. Bias-tees 
are used to split dc bias and the bias-dependent high- 
frequency noise signal. The noise signal is first amplified 
by a low-noise amplifier (LNA) with a noise temperature 
of T^ = 3.5 K in matching conditions, thermalized at 
the same temperature as the sample, then by a series of 
room-temperature amplifiers, and finally collected by a 
zero-bias Schottky diode with band-pass filtering of / = 
600-850 MHz. The noise power measured at the output 
of the LNA is a mixture of thermal noise of the LNA 
and shot noise from the sample. It can be defined in 
terms of the refiected signal |r| (see Fig. 1(b)). Here 
r = fef^ is the signal refiection coefficient when the 
noise generator, i.e. the measured sample with a re- 
sistance i?, does not match to cold amplifier with the 
impedance Zq = 50 Q.. Then, the measured noise power 

reads: P(/) = Pnoz.e(l- |rp) = ^x 2e/x 4Zo (^^ 

where Pnoise = ^2eIR is the shot noise generated by the 
sample at T = 0. When R ^ Zq^ the coupling factor, 
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can be taken as = 1. In our graphene samples, 
this is not fully true, and we employ the differential resis- 
tance measured at bias voltage Vuas for calculating the 
shot noise of the source from the measured noise power. 
The lock-in detection yields the noise derivative Pd = ^, 
which can be integrated to yield the average Fano factor 
F = p^(/o Pddl)/{2el) where Pcai denotes the sepa- 
rately determined calibration constant. Thus, in terms 
of current noise, our average Fano factor corresponds to 
F = (Sjil) - Sm)/{2el) 0. 

The schematic of a typical sample is shown in Fig. 1 (c) . 
Graphene sheets were mechanically exfoliated using the 
Scotch-tape technique [1] and transferred from the natu- 
ral graphite crystals to the surface of Si02/Si substrate 
(300 nm thick thermally grown Si02 layer). The heavily 
doped Si substrate is used as a back-gate in the mea- 
surements. Single graphene layers were located using a 
3CCD camera in an optical microscope on the basis of the 
RGB green component shift, in a similar manner as was 
recently done in Ref. Il5|. After standard e-beam lithogra- 
phy, a bilayer of Ti (10 nm) / Au (40 nm) was evaporated, 
followed by lift-off. Ti was chosen because it allows the 
formation of highly transparent contacts to graphene. In 
fact, in previous work on graphene superconducting junc- 
tions with Al/Ti contacts [16], a very large probability 
for Andreev reffection was observed indicating an average 



transmission probability at the Ti/graphene interface > 
0.8. This is important because in our two lead configu- 
ration, we are sensitive to total transmission eigenvalues 
of both graphene and graphene/leads interface. 
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FIG. 2: Measurements on sample with W jL = 24: (a) Resis- 
tance R (left axis) and conductivity a (right axis) as a func- 
tion of Vgate- (b) Differential resistance dV/dl versus bias 
voltage Vhias at the Dirac point (red curve) and at high den- 
sity (blue curve), (c) Current noise per unit bandwidth Si 
as a function of bias at the Dirac point, at T = 8.5 K, fitted 
(red curve) using Khlus formula (5 = 0.318). Note that the 
low bias data are perfectly fitted as well as the high bias (d) 
Mapping of the average Fano factor F as a function of gate 
voltage Vgate and bias voltage Vhias at T = 8.5 K. 

All the data was measured in a helium dewar, in which 
samples were in a He-gas atmosphere of 1 bar. The resis- 
tance of the samples was measured using standard low- 
frequency ac lock-in technique with an excitation am- 
plitude of 0.3 mV (- 3 K) at uj/2ti = 63.5 Hz, in the 
temperature range of 4.2-30 K. As shown in Fig. 2(a), 
the sample resistance as a function of the charge carrier 
density shows a maximum at the Dirac point for our large 
aspect ratio sample (W/L = 24 with L ~ 200 nm). Our 
graphene samples are non-intentionally p-doped, proba- 
bly due to oxygen gas adsorption [4]. As a consequence, 
the Dirac point is shifted to a positive Vgate value. Never- 
theless, from our measured conductivity values, it seems 
that adsorbed gas on a graphene sheet does not affect 
dramatically the transport properties of our samples. For 
our W/L = 24 sample, we obtain a minimum conductiv- 
ity around ^ which is the one calculated for large aspect 
ratio [5'| and measured in recent experiments ^]. Note 



that the resistance of graphene is nearly independent of 
the bias voltage Vbias (see Fig. 2(b)). 

Fig. 2(c) displays the power spectral density of cur- 
rent noise as a function of the V^ias at T = 8.5 K for 
our W/L = 24 sample, at the Dirac point. The first way 
to extract the Fano factor ^is to use the formula origi- 
nally introduced by Khlus [l7| which describes the cross- 
over from thermal to shot noise when eV ~ ksT [l8|, 

'^^ ~ '^ (r\-z )'^ ^^^^ ( 2k bT ) • Here Zq is the impedance 
of the cold LNA (50 Q). Since the resistance of our 
graphene samples is bias independent, we may fit the 
Khlus formula to our data using only 5^ as a fitting pa- 
rameter at fixed temperature T. From our nearly perfect 
fit, we extract that ^ = 0.318 at T = 9 K . To check 
the accuracy of our measurements, we also used the tun- 
nel junction calibration technique to extract the average 
Fano factor F (8|, (jj. We found F = 0.338 at the Dirac 
point (at Vbias = 40 mV). 

The values found for ^ and F as well as for the mini- 
mum conductivity are very close to the expected theoret- 
ical value of 1/3 and ^ respectively at the Dirac point 
for perfect graphene strip with large W/L [5]. However, 
the Fano factor in this case is also the one expected for 
a diffusive mesoscopic system. In order to confirm that 
charge carriers in our sample do not undergo any scat- 
tering, we swept Vgate to move the Fermi level across the 
Dirac point to high carrier density. A diffusive system 
should not display any gate dependence (as measured by 
DiCarlo et al. in [6| and calculated for long range disor- 
der in [19]). In Fig. 2(d), we can see a mapping of F, 
calculated by integrating the differential Fano factor as 
described in [14], as a function of the Vbias and V„afe • It 
shows a clear dependence of F on gate voltage [20] (i.e. 
the charge carrier density), with a clear drop (about a 
factor of 2) of the Fano factor at large carrier density 
confirming that our results are in good agreement with 
the evanescent state theory [5]. 

It is important to note that the Fano factor is barely af- 
fected by temperature indicating that both the length L 
and the width W are smaller than the electron-phonon 
inelastic scattering length Le-p^: if this condition was 
not fulfilled, the Fano factor would decrease, approxima- 
tively as inversely proportional to N = ^^^xrw-^) ' though 
the actual form would be model dependent [9]. Since our 
shot noise measurements do not depend on temperature 
(between 4 and 30 K) and with our contacts being highly 
transparent, the presence of inelastic scattering mecha- 
nism in the graphene sample and at its interfaces with the 
leads can be ruled out. Note that bad contacts can only 
increase the Fano factor toward the limit of two symmet- 
rical tunneling barriers in series which gives a Fano factor 
of 1/2. This is not the case of our samples in which the 
Fano factor has never been measured higher than 1/3. 

In Fig. 3, we have plotted the F as a function of SV = 
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FIG. 3: F extracted at Vbias = 40 mV for three different 
samples, all having W/L > 3, as a function of SV = Vgate — 
Voirac- For the two unintentionally highly ;?-doped samples 
(orange and green dots), the Dirac point was estimated via 
extrapolation of the minimum conductivity at ^. 



3 (3, 10 and 24). All samples were p-doped, the Dirac 
points being at positive Vgate , but only for one of these 
three samples we could reach the Dirac point (sample 
with W/L = 24). The gate voltages corresponding to 
the Dirac point for the two other samples were estimated 
from their conductivity curves. We observe that despite 
the high doping level of the samples, the Fano factor 
seems to behave universally and tends to zero at very 
high density, i.e. the value for a ballistic system. 

We found that the value for the Fano factor for high 
W/L is about 1/3 to a good accuracy (± 5 %), indicat- 
ing that transport occurs via evanescent modes without 
any observable resonant scattering at the Dirac point. 
It has been theoretically demonstrated that the presence 
of disorder-induced scattering in graphene decreases the 
Fano factor while disorder, count erintuitively, enhances 
conductivity 2l|, |22|. Such behavior may be understood 
as consequences of the absence of intervalley scattering 
[23j and the chirality conservation [2^]. Recently, a model 
calculation near the Dirac point was performed [22] for 
conductivity and shot noise taking into account smooth 
disorder, namely puddles [25|, l26| generated by charged 
impurities. According to this model, smooth disorder at 
length scales <C F, VK lowers the Fano factor at the Dirac 
point, down to 0.243 for one-dimensional disorder and to 
0.295 for the two-dimensional case [22|. This value falls 
clearly below our experimental result ^ = 0.33 ± 0.02 
indicating that transport occurs via evanescent modes 
without any observable scattering at the Dirac point. De- 
spite the probable presence of some disorder in our sys- 
tem, the transport regime can be considered as ballistic 
on our sample length scale. We note that weak disorder 
may induce anomalously large conductance fiuctuations 
J27J. However, we do not observe any oscillations in the 
conductivity or in the Fano factor by increasing density. 

We also measured one sample with a much smaller as- 
pect ratio W/L = 2, having L = 500 nm. The data in 
Fig. 4(a) shows that the minimum conductivity in this 
sample is ~ 6^ ;» ^. The resistance R of the sam- 
ple can be considered as constant as a function of the 
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FIG. 4: Measurements on sample with W/L = 2: (a) i? (left 
axis) and a (right axis) as a function of Vgate- (b) dV/dl 
versus Vbias at the Dirac point (red curve) and at high density 
(blue curve), (c) Current noise per unit bandwidth Si as a 
function of Vbias at the Dirac point, at T = 5 K, fitted (red 
curve) using Khlus formula (5 = 0.196). (d) Mapping of the 
F as a function of Vgate and Vbias at T = 5 K. 



bias (see Fig. 4(b)) In Fig. 4(c), the integrated spectral 
density of current noise as a function of Vbias is fitted 
using Khlus formula. The fit is good at low bias, while 
deviations occur at large Vbias, indicating a decrease of 
the Fano factor, presumably due to elect ron-phonon cou- 
pling [9]. We find that the Fano factor is reduced and 
reaches F = 0.19 at the Dirac point, which is in agree- 
ment with the results of Ref. 5 calculated for the case 
of metallic armchair edge. F as a function of Vbias and 
Vgate is displayed in Fig. 4(d). Note that the determi- 
nation of ^ in Fig. 4(c) yields almost the same result: 
5^ = F = 0.19 at the Dirac point. We observe that F is 
lower than expected from the theory with disorder in Ref. 
2 21 which strengthens the argument that transport at the 
Dirac point can occur without substantial scattering. 

To conclude, we have measured a gate dependent shot 
noise in short graphene strips with large and small W/L. 
At the Dirac point, we observed that for large W/L both 
minimum conductivity and Fano factor reach universal 
values of ^ and 1/3 respectively. For W/L smaller than 
3, the Fano factor is lowered and the minimum conduc- 
tivity increases. These findings are well explained by the 
evanescent wave theory describing transport at the Dirac 
point in perfect graphene. 
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